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ABSTRACT 


Context. Giant radio galaxies (GRGs, or colloquially ‘giants’) are the Universe's largest structures generated by individual galaxies. 
They comprise synchrotron-radiating AGN ejecta and attain cosmological (Mpc-scale) lengths. However, the main mechanisms that 
drive their exceptional growth remain poorly understood. 

Aims. To deduce the main mechanisms that drive a phenomenon, it is usually instructive to study extreme examples. If there exist 
host galaxy characteristics that are an important cause for GRG growth, then the hosts of the largest GRGs are likely to possess 
them. Similarly, if there exist particular large-scale environments that are highly conducive to GRG growth, then the largest GRGs 
are likely to reside in them. For these reasons, we aim to perform a case study of the largest GRG available. 

Methods. We reprocessed the LOFAR Two-metre Sky Survey (LoTSS) DR2 by subtracting compact sources and performing multi- 
scale CLEAN deconvolution at 60" and 90" resolution. The resulting images constitute the most sensitive survey yet for radio 
galaxy lobes, whose diffuse nature and steep synchrotron spectra have allowed them to evade previous detection attempts at higher 
resolution and shorter wavelengths. We visually searched these images for GRGs. 

Results. We discover Alcyoneus, a low-excitation radio galaxy with a projected proper length /, = 4.99 + 0.04 Mpc. Its jets and 
lobes are all four detected at very high significance, and the SDSS-based identification of the host, at spectroscopic redshift Zspec = 
0.24674 + 6-1075, is unambiguous. The total luminosity density at v = 144 MHz is L, = 8+1-107° W Hz !, which is below-average, 
though near-median (percentile 45 + 3%), for GRGs. The host is an elliptical galaxy with a stellar mass M, = 2.4 x 0.4- 10!! Mo and 
a supermassive black hole mass M, = 4+2 -108 Mo, both of which tend towards the lower end of their respective GRG distributions 
(percentiles 25+9% and 23+11%). The host resides in a filament of the Cosmic Web. Through a new Bayesian model for radio galaxy 
lobes in three dimensions, we estimate the pressures in the Mpc?-scale northern and southern lobe to be Pmin,ı = 4.8 +0.3 - 107!6 Pa 
and Pmin2 = 4.9+0.6- 10716 Pa, respectively. The corresponding magnetic field strengths are Bin) = 46+1 pT and Bmin2 = 463 pT. 
Conclusions. We have discovered what is in projection the largest known structure made by a single galaxy — a GRG with a projected 
proper length /, = 4.99 + 0.04 Mpc. The true proper length is at least Imin = 5.04 + 0.05 Mpc. Beyond geometry, Alcyoneus and 
its host are suspiciously ordinary: the total low-frequency luminosity density, stellar mass and supermassive black hole mass are 
all lower than, though similar to, those of the medial GRG. Thus, very massive galaxies or central black holes are not necessary 
to grow large giants, and, if the observed state is representative of the source over its lifetime, neither is high radio power. A low- 
density environment remains a possible explanation. The source resides in a filament of the Cosmic Web, with which it might have 
significant thermodynamic interaction. The pressures in the lobes are the lowest hitherto found, and Alcyoneus therefore represents 
the most promising radio galaxy yet to probe the warm-hot intergalactic medium. 


Key words. galaxies: active — galaxies: individual: Alcyoneus - galaxies: jets - intergalactic medium - radio continuum: galaxies 


1. Introduction 


Most galactic bulges hold a supermassive (ie. M, > 10° Mo) 


Kerr black hole (e.g. [Soltan[1982] that grows by accreting gas, 
dust and stars from its surroundings (Kormendy & Ho|2013). 


The black hole ejects a fraction of its accretion disk plasma 
from the host galaxy along two collimated, magnetised jets 


that are aligned with its rotation axis (e.g. Blandford & Rees 


* [ndear memory of Pallas. If your name hadn't been this popular with 
asteroid discoverers, you'd now be the giants’ giant — once again looking 
down at the sprawling ants below. 


[1974]. The relativistic electrons contained herein experience 
Lorentz force and generate, through spiral motion, synchrotron 
radiation that is observed by radio telescopes. The two jets 
either fade gradually or end in hotspots at the end of diffuse 
lobes, and ultimately enrich the intergalactic medium (IGM) 
with cosmic rays and magnetic fields. The full luminous 
structure is referred to as a radio galaxy (RG). Members of a 
rare RG subpopulation attain megaparsec-scale proper (and 


thus comoving) lengths (e.g. |Willis et al. 
et al. Ishwara-Chandra & Saikia||1999| Jamrozy et al. 


2008| |Machalski Dabhade et al. 
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Fig. 1: Joint radio-infrared view of Alcyoneus, a radio galaxy with a projected proper length of 5.0 Mpc. We show a 
2048" x 2048" solid angle centred around right ascension 123.590372° and declination 52.402795°. We superimpose LOFAR 
Two-metre Sky Survey (LoTSS) DR2 images at 144 MHz of two different resolutions (6” for the core and jets, and 60” for the 
lobes) (orange), with the Wide-field Infrared Survey Explorer (WISE) image at 3.4 um (blue). To highlight the radio emission, the 


infrared emission has been blurred to 0.5’ resolution. 


2020a). The giant radio galaxy (GRG, or colloquially ‘giant’) 
g galaxy quially g 


definition accommodates our limited ability to infer an RG's 
true proper length from observations: an RG is called a GRG 
if and only if its proper length projected onto the plane of the 
sky exceeds some threshold /, Ggc, usually chosen to be 0.7 or 
1 Mpc. Because the conversion between angular length and 
projected proper length depends on cosmological parameters, 
which remain uncertain, it is not always clear whether a given 
observed RG satisfies the GRG definition. 

Currently, there are about a thousand GRGs known, the major- 
ity of which have been found in the Northern Sky. About one 
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hundred exceed 2 Mpc and ten exceed 3 Mpc; at 4.9 Mpc, the 
literature's projectively longest is J1420-0545 
[2008]. As such, GRGs — and the rest of the megaparsec-scale 
RGs — are the largest single-galaxy-induced phenomena in the 
Universe. It is a key open question what physical mechanisms 
lead some RGs to extend for ~10* times their host galaxy 
diameter. To determine whether there exist particular host 
galaxy characteristics or large-scale environments that are 
essential for GRG growth, it is instructive to analyse the largest 
GRGs, since in these systems it is most likely that all major 
favourable growth factors are present. We thus aim to perform 


Martijn S.S.L. Oei et al.: The discovery of a radio galaxy of at least 5 Mpc 


a case study of the largest GRG available. 


As demonstrated by |Dabhade et al.| (2020b)’s record sample 
e Low- 


of 225 discoveries, t requency Array (LOFAR) 
is among the most attractive contempo- 
rary instruments for finding new GRGs. This Pan-European 
radio interferometer features a unique combination of short 
baselines to provide sensitivity to large-scale emission, and 
long baselines to mitigate source confusion[!] These qualities 
are indispensable for observational studies of GRGs, which 
require identifying both extended lobes and compact cores and 
jets. Additionally, the metre wavelengths at which the LOFAR 
operates allow it to detect steep-spectrum lobes far away from 
host galaxies. Such lobes reveal the full extent of GRGs, but 
are missed by decimetre observatories. Thus, in Section [2] we 
describe a reprocessing of the LOFAR Two-metre Sky Survey 
(LoTSS) Data Release 2 (DR2) aimed at revealing hitherto 
unknown RG lobes — among other goals. An overview of the 
reprocessed images, which cover thousands of square degrees, 
and statistics of the lengths and environments of the GRGs 
they have revealed, are subjects of future publications. For 
now, these images allow us to discover Alcyoneuf] a 5 Mpc 
GRG, whose properties we determine and discuss in Section|3| 
Figure |1| provides a multi-wavelength, multi-resolution view 
of this giant. ose ius our concluding remarks. 

We assume a concordance inflationary ACDM model with 
parameters Ji from[Planck Collaboration et al.| (2020); i.e. Mt = 
(h = 0.6766, Og = 0.0490, Quo = 0.3111, Q4 o = 0.6889), 
where Ho := h- 100 kms ! Mpc™!. We define the spectral 
index a such that it relates to flux density F, at frequency v 
as F, œ v”. Regarding terminology, we strictly distinguish 
between a radio galaxy, a radio-bright structure of relativistic 
particles and magnetic fields (consisting of a core, jets, hotspots 
and lobes), and the host galaxy that generates it. 


2. Data and methods 
The LoTSS, conducted by the LOFAR High-band Antennae 


(HBA), is a 120-168 MHz interferometric survey (Shimwell et al. 
2017||2019 with the ultimate aim to image the full 


Northern Sky at resolutions of 6", 20", 60" and 90". Its cen- 
tral frequency v; = 144 MHz. The latest data release — the 
LoTSS DR2 — covers 27% of the North- 
ern Sky, split over two regions of 4178 deg? and 1457 deg?; the 
largest of these contains the Sloan Digital Sky Survey (SDSS) 
DR7 area. By default, the LoTSS DR2 pro- 
vides imagery at the 6" and 20" resolutions. We show these 
standard products in Figure [2] for the same sky region as in 
Figure |1| In terms of total source counts, the LoTSS DR2 is 
the largest radio survey carried out thus far: its catalogue con- 
tains 4.4 - 10° sources, most of which are considered com- 
pact. By contrast, the 60" and 90" imagery, which we dis- 


cuss in more detail in (in prep}. is intended to re- 


veal extended structures in the low-frequency radio sky, such 


! Source confusion is an instrumental limitation that arises when the 
resolution of an image is low compared to the sky density of statisti- 
cally significant sources. It causes angularly adjacent, but physically 


unrelated sources to blend together, making it hard or even impossible 
to distinguish them (e.g.|Condon et al.|2012). 
? Alcyoneus was the son of Ouranos, the Greek primordial god of the 


sky. According to Ps.-Apollodorus, he was also one of the greatest of 
the Gigantes (Giants), and a challenger to Heracles during the Gigan- 
tomachy — the battle between the Giants and the Olympian gods for 
supremacy over the Cosmos. The poet Pindar described him as ‘huge 
as a mountain’, fighting by hurling rocks at his foes. 
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Fig. 2: Aleyoneus’ lobes are easily overlooked in the LoTSS 
DR2 at its standard resolutions. We show images at cen- 
tral frequency v; = 144 MHz and resolutions rwym = 6" 
(top) and 6gwum = 20" (bottom), centred around host galaxy 
J081421.68--522410.0. 


as giant radio galaxies, supernova remnants in the Milky Way, 
radio halos and shocks in galaxy clusters, and — potentially 
— accretion shocks or volume-filling emission from filaments 
of the Cosmic Web. To avoid the source confusion limit at 
these resolutions, ee eres cL we used 
DDFacet (Tasse et al.|2018) to predict visibilities corresponding 
to the 20" LoTSS DR2 sky model and subtracted these from the 
data, before imaging at 60" and 90" with WSClean IDG 
. We used -0.5 Briggs 
weighting and multiscale CLEAN (Offringa & Smirnov]|2017}, 
with -multiscale-scales 0,4,8,16,32,64. Importantly, we 
did not impose an inner (u, v)-cut. We imaged each pointing sep- 
arately, then combined the partially overlapping images into a 
mosaic by calculating, for each direction, a beam-weighted av- 
erage. 

Finally, we visually searched the LoTSS DR2 for GRGs, primarily 
at 6" and 60" using the Hierarchical Progressive Survey (HiPS) 


system in Aladin Desktop 11.0 (Bonnarel et al.|2000). 
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3. Results and discussion 
3.1. Radio morphology and interpretation 


During our LoTSS DR2 search, we identified a three-component 
radio structure of total angular length ø = 20.8’, visible at all 
(6", 20", 60” and 90") resolutions. Figure|2|provides a sense of 
our data quality; it shows that the outer components are barely 
discernible in the LoTSS DR2 at its standard 6” and 20” reso- 
lutions. Meanwhile, Figure [1| shows the outer components at 
60", and the top panel of Figure[|9|shows them at 90"; at these 
resolutions, they lie firmly above the noise. Compared with the 
outer structures, the central structure is bright and elongated, 
with a 155" major axis and a 20" minor axis. The outer struc- 
tures lie along the major axis at similar distances from the cen- 
tral structure, are diffuse and amorphous, and feature specific 
intensity maxima along this axis. 

In the arcminute-scale vicinity of the outer structures, the DESI 


Legacy Imaging Surveys (Dey et al.|2019) DR9 does not reveal 


galaxy overdensities or low-redshift spiral galaxies, the ROSAT 


All-sky Survey (RASS) (Voges et al.|1999) does not show X-ray 


brightness above the noise, and there is no Planck Sunyaev- 
Zeldovich catalogue 2 (PSZ2) 
source nearby. The outer structures therefore cannot be super- 
nova remnants in low-redshift spiral galaxies or radio relics 
and radio halos in galaxy clusters. Instead, the outer structures 
presumably represent radio galaxy emission. The radio-optical 
overlays in Figure|3|s top and bottom panel show that it is im- 
probable that each outer structure is a radio galaxy of its own, 
given the lack of significant 6" radio emission (solid light green 
contours) around host galaxy candidates suggested by the mor- 
phology of the 60" radio emission (translucent white contours). 
For these reasons, we interpret the central (jet-like) structure 
and the outer (lobe-like) structures as components of the same 
radio galaxy. 

Subsequent analysis — presented below — demonstrates that 
this radio galaxy is the largest hitherto discovered, with a pro- 
jected proper length of 5.0 Mpc. We dub this GRG Alcyoneus. 


3.2. Host galaxy identification 


Based on the middle panel of Figure[3|and an SDSS DR12 
et al.[2015) spectrum, we identify a source at a J2000 right as- 
cension of 123.590377^, a declination of 52.402795° and a spec- 
troscopic redshift of Zspec = 0.24674 + 6- 107? as Alcyoneus' 
host. Like most GRG hosts, this source, with SDSS DR12 name 
J081421.68+522410.0, is an elliptical galaxyf] without a quasar. 
From optical contours, we find that the galaxy’s minor axis 
makes a ~20° angle with Alcyoneus’ jet axis. 

In Figure we further explore the connection between 
J081421.68+522410.0 and Alcyoneus’ radio core and jets. From 
top to bottom, we show the LoTSS DR2 at 6”, the Very Large 


Array Sky Survey (VLASS) at 2.2", and the 
an 


cy et al. 
Panoramic Survey Telescope and Rapid Response System (Pan- 
STARRS) DR1 i-band. Two facts confirm 
that the host identification is highly certain. First, for both the 
LoTSS DR2 at 6” and the VLASS at 2.2”, the angular separa- 
tion between J081421.68+522410.0 and the arc connecting Al- 
cyoneus' two innermost jet features is subarcsecond. Moreover, 
the alleged host galaxy is the brightest Pan-STARRS DRI i-band 


? Based on the SDSS morphology,|Kuminski & Shamir calculate 


a probability of 89% that the galaxy is an elliptical. 
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Fig. 3: Joint radio-optical views show that Figure [1]s outer 
structures are best interpreted as a pair of radio galaxy 
lobes fed by central jets. On top of DESI Legacy Imaging 
Surveys DR9 (g, r, z)-imagery, we show the LoTSS DR2 at var- 
ious resolutions through contours at multiples of c, where c 
is the image noise at the relevant resolution. The top and bot- 
tom panel show translucent white 60" contours at 3,5, 7,9, 110 
and solid light green 6" contours at 4, 7, 10, 20, 400. The central 
panel shows translucent white 6” contours at 5, 10, 20, 40, 800. 
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Fig. 4: The SDSS DR12 source J081421.68+522410.0 is Al- 
cyoneus' host galaxy. The panels cover a 2.5' x 2.5' region 
around J081421.68+522410.0, an elliptical galaxy with spectro- 
scopic redshift Zspec = 0.24674 + 6- 1075. From top to bot- 
tom, we show the LoTSS DR2 6", the VLASS 2.2”, and the Pan- 
STARRS DRI i-band — relative to the peak specific intensity of 
J081421.68+522410.0 — with LoTSS contours (white) as in Fig- 
ure[3|and a VLASS contour (gold) at 5c. 


source within a radius of 45" of the central VLASS image com- 
ponent. 


3.3. Radiative- or jet-mode active galactic nucleus 


Current understanding (e.g. Heckman & Best||2014) suggests 
that the population of active galactic nuclei (AGN) exhibits a 


dichotomy: AGN seem to be either radiative-mode AGN, which 
generate high-excitation radio galaxies (HERGs), or jet-mode 
AGN, which generate low-excitation radio galaxies (LERGs). Is 
Alcyoneus a HERG or a LERG? The SDSS spectrum of the host 
features very weak emission lines; indeed, the star formation 
rate (SFR) is just 1.6 - 107? Mo yr^! (Chang et al.||2015). Fol- 
lowing the classification rule of|Best & Heckman] (2012); 
let al. (2014); [Pracy et al) (2016); 2018) based on 


the strength and equivalent width of the OIII 5007 À line, we 
conclude that Alcyoneus is a LERG. Moreover, the WISE pho- 


tometry (Cutri & etal.|2012) at 11.6 um and 22.1 um is below the 


instrumental detection limit. Following the classification rule of 


based on the 22.1 um luminosity density, 
we affirm that Alc 


yoneus is a LERG. Through automated classi- 
fication, [Best & Heckman]|(2012) came to the same conclusion. 
Being a jet-mode AGN, the supermassive black hole (SMBH) in 
the centre of Alcyoneus' host galaxy presumably accretes at an 
efficiency below 146 of the Eddington limit, and is fueled mainly 
by slowly cooling hot gas. 


3.4. Projected proper length 


We calculate Alcyoneus' projected proper length /, through 
its angular length $ and spectroscopic redshift Zspec. We for- 
mally determine ¢ = 20.8’ + 0.15’ from the compact-source- 
subtracted 90” image (top panel of Figure |9) by selecting the 
largest great-circle distance between all possible pairs of pix- 
els with a specific intensity higher than three sigma-clipped 
standard deviations above the sigma-clipped median. We find 
lp = 4.99 + 0.04 Mpc; this makes Alcyoneus the projectively 
largest radio galaxy known. For methodology details, and for a 
probabilistic comparison between the projected proper lengths 
of Alcyoneus and J1420-0545, see Appendix[A] 


3.5. Radio luminosity densities and kinetic jet powers 


From the LoTSS DR2 6" image (top panel of Figure|4), we mea- 
sure that two northern jet local maxima occur at angular dis- 
tances of 9.2 + 0.2” and 23.7 + 0.2” from the host, or at pro- 
jected proper distances of 36.8 + 0.8 kpc and 94.8 + 0.8 kpc. 
Two southern jet local maxima occur at angular distances of 
8.8 + 0.2” and 62.5 + 0.2" from the host, or at projected proper 
distances of 35.2 + 0.8 kpc and 249.9 + 0.8 kpc. 

At the central observing frequency of v, = 144 MHz, the north- 
ern jet has a flux density F, = 193+20 mJy, the southern jet has 
F, = 110x12 mJy, whilst the northern lobe has F, = 63-7 mJy 
and the southern lobe has F, = 44 + 5 mJy. To minimise con- 
tamination from fore- and background galaxies, we determined 
the lobe flux densities from the compact-source-subtracted 90" 
image. The flux density uncertainties are dominated by the 1096 
flux scale uncertainty inherent to the LoTSS DR2 
. The host galaxy flux density is relatively weak, 
and the corresponding emission has, at v; = 144 MHz and 
6" resolution, no clear angular separation from the inner jets' 
emission; we have therefore not determined it. 

Due to cosmological redshifting, the conversion between flux 
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Fig. 5: The LoTSS-VLASS spectral index map reveals Alcy- 
oneus' flat-spectrum core and steeper-spectrum jets. We 
show all directions where both the LoTSS and VLASS image 
have at least 50 significance. In black, we overlay the same 
LoTSS contours as in Figures 3|and|4 4 The core spectral index 


is a = —0.25 + 0.1 and the combined inner jet spectral index is 
a = —0.65 x 0.1. 


density and luminosity density depends on the spectral indices 
a of Alcyoneus’ luminous components. We estimate the spec- 
tral indices of the core and jets from the LoTSS DR2 6" and 
VLASS 2.2" images. After convolving the VLASS image with 
a Gaussian to the common resolution of 6", we calculate the 
mean spectral index between LoTSS' v, = 144 MHz and VLASS’ 
y, = 2.99 GHz. Using only directions for which both images 
have a significance of at least 5c", we deduce a core spectral 
index a = —0.25 + 0.1 and a combined inner jet spectral index 
a = —0.65 x0.1. The spectral index uncertainties are dominated 
by the LoTSS DR2 and VLASS flux scale uncertainties. We show 
the full spectral index map in Figure[5| We have not determined 
the spectral index of the lobes, as they are only detected in the 
LoTSS imagery. 

The luminosity densities of the northern and southern jet at 
rest-frame frequency v = 144 MHz are L, = (3.6 + 0.4). 
10? W Hz ! and L, = (2.0 + 0.2) - 10? W Hz |, respectively. 


Following|Dabhade et al.|(2020a), we estimate the kinetic power 


of the jets from their luminosity densities and the e of the 
simulation-based analytical model of[Hardcastle|(2018) 2018). We find 
Oye. = 1.2 € 0.1- 106 W and Qiu2 = 66 + A 10? W, 
so that the total kinetic jet power is Qjets = Qjet1 + Qjet2 = 
1.9 + 0.2 - 1076 W. Interestingly, this total kinetic jet power is 
lower than the average Qjets = 3.7 - 10?6 W, and close to the me- 
dian Qia, = 2.2- 10°° W, for low-excitation giant radio galaxies 


(LEGRGs) in the redshift range 0.18 « z < 0.43 (Dabhade et al. (Dabhade et al.] 
Because the lobe spectral indices are unknown, we present lu- 
minosity densities for several possible values of «œ in Table[1]] 
(Because of electron ageing, œ will decrease further away from 
the core.) 


* 'The inferred luminosity densities have a cosmology-dependence; 
our results are ~6% higher than for modern high-H, cosmologies. 
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Fig. 6: Alcyoneus has a low-frequency luminosity density 
typical for GRGs. We explore the relation between GRG pro- 
jected proper length /, and total luminosity density L, at rest- 
frame frequency v = 144 MHz. Total luminosity densities in- 
clude contributions from all available radio galaxy components 
(i.e. the core, jets, hotspots and lobes). Literature GRGs are from 
(2020b), and are marked with grey disks, while 
Alcyoneus is marked with a green star. Translucent ellipses in- 
dicate -1 to +1 standard deviation uncertainties. Alcyoneus has 
a typical luminosity density (percentile 45 + 3%). 


Table 1: Luminosity densities L, (in 1074 W Hz!) of Alcyoneus' 
lobes for three potential spectral indices œ at rest-frame fre- 


quency v = 144 MHz, assuming a|Planck Collaboration et al. 
cosmology. 


a = —0.8 
12+1 
8.3 + 0.8 


y--12 a=-1.6 
1321 14+1 
9.0 x 0.9 9.04] 


Northern lobe 


Southern lobe 


Assuming œ = -—1.2, Alcyoneus total luminosity density at 
v = 144 MHz is L, = 7.8 + 0.8 - 10? WHz !. In Figure 
compare this estimate to other GRGs' total luminosity density at 


the same E as found i. (2020b) through 
the LoTSS DR1 (Shimwell et al.[2019). Interestingly, Alcyoneus 
is not particularly luminous: T as a low-frequency luminosity 
density typical for the currently known GRG population (per- 
centile 45 + 3%). 


3.6. True proper length: relativistic beaming 


Following |Hardcastle et al.| (1998a), we simultaneously con- 


strain Alcyoneus' jet speed u and inclination angle 0 from the 
jets’ flux density asymmetry: the northern-to-southern jet flux 
density ratio J = 1.78 + 0.3] We assume that the jets prop- 
agate with identical speeds u in exactly opposing directions 
(making angles with the line-of-sight 0 and 0 + 180°), and have 
statistically identical relativistic electron populations, so that 
they have a common synchrotron spectral index a. Using a = 
—0.65 + 0.1 as before, and 


u Ja — 1 
B; Bose = ——— 
C JZ +i 


(1) 


A i 
- vL 


? Because J is obtained through division of two independent normal 
random variables (RVs) with non-zero mean, J is an RV with an un- 
correlated noncentral normal ratio distribution. 
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we find B cos 0 = 0.106 + 0.03. Because cos 0 x 1, B is bounded 
from below by Bmin = 0.106 + 0.03. 

From detailed modelling of ten Fanaroff-Riley (FR) I radio 
aN (which have jet luminosities comparable to Alcy- 


oneus’ ), {Laing & Bridle|(2014) 2014) deduced that initial jet speeds are 
roughly B = 0.8, which decrease until roughly 0.6 ro, with ro be- 
ing the recollimation distance. Most of UTE 
ten recollimation distances are between 5 and 15 kpc, with the 
largest being that of NGC 315: rọ = 35 kpc. Because the lo- 
cal specific intensity maxima in Alcyoneus' jets closest to the 
host occur at projected proper distances of 36.8 + 0.8 kpc and 
35.2 0.8 kpc, the true proper distances must be even larger. We 
conclude that the observed jet emission presumably comes from 
a region further from the host than ro, so that the initial stage 
of jet deceleration — in which the jet speed is typically reduced 
by several tens of percents of c — must already be completed. 
Thus, Bmax = 0.8 is a safe upper bound. 

Taking Bmax = 0.8, 0 is bounded from above by max = 82.4 + 2° 
(0 € [0,90°]), or bounded from below by 180° — Omax = 97.6 + 
2° (0 e [90, 180*f If we model Alcyoneus' geometry as a 
line segment, and assume no jet reorientation, Alcyoneus' true 
proper length / and projected proper length /, relate as 


M NM (2) 


^ SinÓg 

We bound / from below: min = 5.04 + 0.05 Mpc. A triangu- 
lar prior on B between Bmin and Bmax with the mode at Bmax 
induces a skewed prior on /; the 90% credible interval is / € 
[5.0 Mpc, 5.5 Mpc], with the mean and median being 5.2 Mpc 
and 5.1 Mpc, respectively. A flat prior on B between fj, and 
Bmax also induces a skewed prior on /; the 90% credible inter- 
val is / € [5.0 Mpc, 7.1 Mpc], with the mean and median being 
5.6 Mpc and 5.1 Mpc, respectively. The median of / seems par- 
ticularly well determined, as it is insensitive to variations of the 
prior on f. 

In Appendix[B| we explore the inclination angle conditions un- 
der which Alcyoneus has the largest true proper length of all 
known (» 4 Mpc) GRGs. 


3.7. Stellar and supermassive black hole mass 


Does a galaxy or its central black hole need to be massive in 
order to generate a GRG? 

Alcyoneus' host has a stellar mass M, = 2.4 € 0.4 - 10!! M 
(Chang et al.|2015). We test whether or not this is a typical stel- 
lar mass among the total known GRG population. We assem- 
ble a literature catalogue of 1013 GRGs by merging the com- 
pendium of |Dabhade et al.| (2020a), which is = up to 
April 2020, with the GRGs icon = 


Ishwara-Chandra et aL 2020) [T 
br Rater ey Balhae eral ]po21)] 


Eoi 
Mahato et al.|(2021). We ees sieur masses with hies 


from 1.|(2015] , which 2 based on SDSS and WISE 
photometry, and trom 1.) (2018), which are based on 
GALEX, SDSS and WISE photometry We give precedence to 
the stellar masses by |S 8) when both are avail- 
able. We obtain stellar masses for =e a ously known GRGs. 
The typical stellar mass range is 10!! - 10!? Mo, the median 


$ Taking Bmax = 1 instead, 6 is bounded from above by Omax = 83.9+2° 
(0 € [0,90°]), or bounded from below by 180° — Omax = 96.1 + 2° (0 € 
[90°, 180*]). 
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Fig. 7: Alcyoneus' host has a lower stellar and supermas- 
sive black hole mass than most GRG hosts. We explore 
the relations between GRG projected proper length /, and host 
galaxy stellar mass M, (top panel) or host galaxy supermas- 
sive black hole mass M, (bottom panel). Our methods allow de- 
termining these properties for a small proportion of all litera- 
ture GRGs only. Literature GRGs are marked with grey disks, 
while Alcyoneus is marked with a green star. Translucent el- 
lipses indicate -1 to +1 standard deviation uncertainties. Alcy- 
oneus' host has a fairly typical — though below-average — stel- 
lar mass (percentile 25 x 996) and supermassive black hole mass 
(percentile 23 + 11%). 


M, = 3.5 - 10!! M, and the mean M, = 3.8 - 10!! Mo. Strik- 
ingly, the top panel of Figure[|illustrates that Alcyoneus' host 
has a fairly low (percentile 25 +9%) stellar mass compared with 
the currently known population of GRG hosts. 

For the GRGs in our literature catalogue, we also estimate 
SMBH masses via the M-sigma relation. We collect SDSS DR12 


stellar velocity dispersions with uncertainties (Alam et al.|2015), 
and apply the M-sigma relation of Equation 7 in 
[& Ho| (2013). Alcyoneus’ host has a SMBH mass M, = 3.9 + 
1.7-10° Mo. We obtain SMBH masses for 189 previously known 
GRGs. The typical SMBH mass range is 108 — 10? Mo, the me- 
dian M, = 7.9- 108 Mo and the mean M, = 1.5 - 10? Mo. 
Strikingly, the bottom panel of Figure []illustrates that Alcy- 
oneus' host has a fairly low (percentile 23 + 11%) SMBH mass 
compared with the currently known population of GRG hosts. 

We note that Alcyoneus is the only GRG with l, > 3 Mpc whose 
host's stellar mass is known through(Chang etal 2015) o a 
be estimated 


, and whose host’s SMBH mass can 
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through its SDSS DR12 velocity dispersion. These data allow 
us to state confidently that exceptionally high stellar or SMBH 
masses are not necessary to generate 5-Mpc-scale GRGs. 


3.8. Surrounding large-scale structure 


Several approaches to large-scale structure (LSS) classification, 
such as the T-web scheme (Hahn et al.[2007], partition the mod- 
ern Universe into galaxy clusters, filaments, sheets and voids. In 
this section, we determine Alcyoneus' most likely environment 
type. 
We conduct a tentative quantitative analysis using the SDSS 
DR7 spectroscopic galaxy sample ATETA Does 
Alcyoneus’ host have fewer, about equal or more galactic neigh- 
bours in SDSS DR7 than a randomly drawn galaxy of similar 
r-band luminosity density and redshift? Let r (z) be the comov- 
ing radial distance corresponding to cosmological redshift z. We 
consider a spherical shell with the observer at the centre, inner 
radius max (r(z = Zspec) — ro, 0) and outer radius r(z = Zspec) + ro. 
We approximate Alcyoneus' cosmological redshift with Zgpec 
and choose rọ = 25 Mpc. As all galaxies in the spherical shell 
have a similar distance to the observer (i.e. distances are at most 
2ro different), the SDSS DR7 galaxy number density complete- 
ness must also be similar throughout the spherical shell|’| For 
each enclosed galaxy with an r-band luminosity density be- 
tween | — 6 and 1 + 6 times that of Alcyoneus’ host, we count 
the number of SDSS DR7 galaxies N<pr (R) within a sphere of 
comoving radius R around it — regardless of luminosity den- 
sity, and excluding itself. Alcyoneus' host has an SDSS r-band 
apparent magnitude m, — 18.20; the corresponding luminosity 
density is L, (A, = 623.1 nm) = 3.75 - 10°? W Hz^!. We choose 
6 = 0.25; this yields 9,358 such enclosed galaxies. 
In d we show the distribution of Ner(R) for R = 
5 Mpc and R = 10 Mpc. We verify that the distributions 
are insensitive to reasonable changes in ro and ô. Note that 
there is no SDSS DR7 galaxy within a comoving distance 
of 5 Mpc from Alcyoneus' host. The nearest such galaxy, 
J081323.49+524856.1, occurs at a comoving distance of 7.9 Mpc: 
the nearest ~2, 000 Mpc? of comoving space are free of galactic 
neighbours with L, (A,) > 5.57 -10° W Hz ! Phin the same way 
as in Section[.1] we verify that the DESI Legacy Imaging Sur- 
veys DR9, RASS and PSZ2 do not contain evidence for a galaxy 
cluster in the direction of Alcyoneus' host. The nearest galaxy 
cluster, according to the SDSS-III cluster catalogue of|Wen et al. 
(2012), instead lies 24’ away at right ascension 123.19926’, dec- 
ination 52.72468° and photometric redshift zp, = 0.2488. It has 
an Rogo = 1.1 Mpc and, according to the DESI cluster catalogue 
ofouet al nan atotal mass M = 2.2-10'* Mo. The comoving 
distance between the cluster and Alcyoneus' host is 11 Mpc. All 
in all, we conclude that Alcyoneus does not reside in a galaxy 
cluster. Meanwhile, there are five SDSS DR7 galaxies within a 
comoving distance of 10 Mpc from Alcyoneus' host: this makes 
it implausible that Alcyoneus lies in a void. Finally, one could 
interpret Nr (R) as a proxy for the LSS total matter density 
around a galaxy. For R = 10 Mpc, just 17% of galaxies in the 
shell with a similar luminosity density as Alcyoneus’ host have 
a higher LSS total matter density. Being on the high end of the 


7 For ro = 25 Mpo, this is a good approximation, because the shell is 
cosmologically thin: 279 = 50 Mpc roughly amounts to the length of a 
single Cosmic Web filament. 

8 This is the luminosity density that corresponds to the SDSS r-band 
apparent magnitude completeness limit m, = 17.77 (Strauss et al. 


poa 
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Fig. 8: Like most galaxies of similar r-band luminosity 
density and redshift, Alcyoneus' host has no galactic 
neighbours in SDSS DR7 within 5 Mpc. However, within 
10 Mpc, Alcyoneus’ host has more neighbours than most 
similar galaxies. For all 9,358 SDSS DR7 galaxies with an r- 
band luminosity density between 75% and 125% that of AI- 
cyoneus’ host and a comoving radial distance that differs at 
most 79 = 25 Mpc from Alcyoneus’, we count the number of 
SDSS DR7 galaxies Nep (R) within a sphere of comoving radius 
R = 5 Mpc (top panel) and R = 10 Mpc (bottom panel). The top 
panel indicates that Alcyoneus does not inhabit a galaxy clus- 
ter; the bottom panel indicates that Alcyoneus does not inhabit 
a void. 


density distribution, but lying outside a cluster, Alcyoneus most 
probably inhabits a filament of the Cosmic Web. 


3.9. Proper lobe volumes 


We determine the proper volumes of Alcyoneus' lobes with a 
new Bayesian model. The model describes the lobes through a 
pair of doubly truncated, optically thin cones, each of which 
has a spatially constant and isotropic monochromatic emis- 
sion coefficient (MEC) GEN m T We allow 
the 3D orientations and opening angles of the cones to differ, 
as the lobes may traverse their way through differently pres- 
sured parts of the warm-hot intergalactic medium (WHIM): 
e.g. the medium near the filament axis, and the medium near 
the surrounding voids. By adopting a spatially constant MEC, 
we neglect electron density and magnetic field inhomogeneities 
as well as spectral-ageing gradients; by adopting an isotropic 
MEC, we assume non-relativistic velocities within the lobe so 
that beaming effects are negligible. Numerically, we first gener- 
ate the GRG's 3D MEC field over a cubical voxel grid, and then 
calculate the corresponding model image through projection, 
including expansion-related cosmological effects. Before com- 
parison with the observed image, we convolve the model image 
with a Gaussian kernel to the appropriate resolution. We exploit 
the approximately Gaussian LoTSS DR2 image noise to formu- 
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late the likelihood, and assume a flat prior distribution over the 
parameters. Using a Metropolis-Hastings (MH) Markov chain 
Monte Carlo (MCMC), we sample from the posterior distribu- 
tion 

In the top panel of Figure(9 we show the LoTSS DR2 compact- 
source-subtracted 90” image of Alcyoneus. The central region 
has been excluded from source subtraction, and hence Alcy- 
oneus’ core and jets remain. (However, when we run our MH 
MCMC on this image, we do mask this central region.) In the 
middle panel, we show the highest-likelihood (and thus max- 
imum a posteriori (MAP)) model image before convolution. In 
the bottom panel, we show the same model image convolved 
to 90" resolution, with 2c and 30 contours of the observed im- 
age overlaid. We provide the full parameter set that corresponds 
with this model in Table[C.1] 

The posterior mean, calculated through the MH MCMC sam- 
ples after burn-in, suggests the following geometry. The north- 
ern lobe has an opening angle y, = 10 + 1°, and the cone 
truncates at an inner distance dj; = 2.6 + 0.2 Mpc and at an 
outer distance do, = 4.0 + 0.2 Mpc from the host galaxy. The 
southern lobe has a larger opening angle y2 = 26 + 2°, but its 
cone truncates at smaller distances of di2 = 1.5 + 0.1 Mpc and 
do» = 2.0 x 0.1 Mpc from the host galaxy. These parameters fix 
the proper volumes of Alcyoneus' northern and southern lobes. 
We find Vj = 1.5 + 0.2 Mpc? and V; = 1.0+0.2 Mpc’, respec- 
tively (see Equation|C.15} 

How are the lobes oriented? Figure [1] provides a visual hint 
that the lobes are subtly non-coaxial. The posterior indicates 
that the position angles of the northern and southern lobes are 
gı = 307+1° and y» = 139+2’, respectively. The position angle 
difference is thus Ay = 168+2°: although close to Ay = 180°, we 
can reject coaxiality with high significance. Interestingly, the 
posterior also constrains the angles that the lobe axes make with 
the plane of the sky: |8; — 90°| = 51 + 2° and |@, - 90°| = 18 + 7°. 
Again, the uncertainties imply that the lobes are probably not 
coaxial. We stress that these inclination angle results are tenta- 
tive only. Future model extensions should explore how sensitive 
they are to the assumed lobe geometry (by testing other shapes 
than just truncated cones, such as ellipsoids). 

One way to validate the model is to compare the observed 
lobe flux densities of Section to the predicted lobe flux 
densities. According to the posterior, the MECs of the north- 
ern and southern lobes are j,,; = 17 + 2 Jy deg ? Mpc'! and 
j2 = 18 +3 Jy deg ? Mpc !. Combining MECs and volumes, 
we predict northern and southern lobe flux densities F,(v.) = 
63 + 4 mJy and Fy5(v.) = 45 + 5 mJy (see Equation E 
We find excellent agreement: the relative differences with the 
observed results are 0% and 2%, respectively. 


? For a detailed description of the model parameters, the MH MCMC 
and formulae for derived quantities, see Appendix 

? As a sanity check, we compare our results to those from a less rig- 
orous, though simpler ellipsoid-based method of estimating volumes. 
By fitting ellipses to Figure [9}s top panel image, one obtains a semi- 
minor and semi-major axis; the half-diameter along the ellipsoid’s third 
dimension is assumed to be their mean. This method suggests a north- 
ern lobe volume V; = 1.4 + 0.3 Mpc? and a southern lobe volume 
V; = 1.1 + 0.3 Mpc’. These results agree well with our Bayesian 
model results. (If the half-diameter along the third dimension is instead 
treated as an RV with a uniform distribution between the semi-minor 
axis and the semi-major axis, the estimates remain the same.) 
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Fig. 9: Alcyoneus' lobe volumes can be estimated by com- 
paring the observed radio image to modelled radio im- 
ages. Top: LoTSS DR2 compact-source-subtracted 90" image of 
Alcyoneus. For scale, we show the stellar Milky Way disk (di- 
ameter: 50 kpc) and a 10 times inflated version; the spiral galaxy 
shape D ME MEN EE Middle: Highest- 
likelihood model image. Bottom: The same model image con- 
volved to 90" resolution, with 207 and 30 contours of the ob- 


served image overlaid. Article number, page 9 of 18 


A&A proofs: manuscript no. aanda 


3.10. Lobe pressures and the local WHIM 


From Alcyoneus’ lobe flux densities and volumes, we can in- 
fer lobe pressures and magnetic field strengths. We calculate 
these through pysyncH | , which uses 


(Hardcastle et al. 
the formulae first proposed by|Myers & Spangler and re- 
examined by |Beck & Krause} (2005). Following the notation of 


Hardcastle et al.|(1998b), we assume that the electron energy 
distribution is a power law in Lorentz factor y with ymin = 10, 
Ymax = 10* and exponent p — —2; we also assume that the 
kinetic energy density of protons is vanishingly small com- 
pared with that of electrons (x = 0), and that the plasma fill- 
ing factor is unity ( = 1). Assuming the minimum-energy 
condition (Burbidge[1956], we find minimum-energy pressures 
Prin.) = 4.8+0.3-10 © Pa and Pyin = 4.9+0.6-10-'° Pa for the 
northern and southern lobes, respectively. The corresponding 
minimum-energy magnetic field strengths are Bmin,ı = 4641 pT 
and Bmin2 = 46 + 3 pT. Assuming the equipartition condi- 
tion SEES pori], we find equipartition pressures Pog = 
4.9 + 0.3 - 10 * Pa and Peq2 = 4.9 + 0.6 - 10 !6 Pa for the 
northern and southern lobes, respectively. The corresponding 
equipartition magnetic field strengths are Beg, = 43 +2 pT and 
Beg = 43 +2 pT. The minimum-energy and equipartition re- 
sults do not differ significantly. 

From pressures and volumes, we estimate the internal energy 
of the lobes E = 3PV. We find Egi; = 6.2 + 0.5 - 10? J, 
Emin2 = 4.3 + 0.6- 10 J, E41 = 63 + 0.5- 10°? J and Enz = 
4.4 + 0.6 - 10° J. Next, we can bound the ages of the lobes from 
below by neglecting synchrotron losses, and assuming that the 
jets have been injecting energy in the lobes continuously at the 
currently observed kinetic jet powers. Using At = E n we 
find Atmin, = 1.7 € 0.2 Gyr, Atmin = 2.1 + 0.4 Gyr, and identi- 
cal results when assuming the equipartition condition. Finally, 
we can obtain a rough estimate of the average expansion speed 
of the radio galaxy during its lifetime u = I,(Ar) |. We find 
u = 2.6 + 0.3- 10? km sl, or about 1% of the speed of light. 


Several other authors 
Subrahmanyan et al.[1996]|Parma et aL[1996] [Mack et aL [998 
Schoenmakers et al.|1998/2000||shwara-Chandra & Saikia|1999 
Tara et al. [2000] [Machalski & Jamrozy[2000] [Machalski et al. 
2001] Saripalli et al[2002][amrozy et al.]2005| Bubrahmanyan 
et al.|2006][2008] Saikia et aL[2006]|Machalski et aL.[2006][2007 
Safouris et al.|2009}|Malarecki et al.]2013 


have estimated the minimum-energy or equipartition 
pressure of the lobes of GRGs embedded in non-cluster envi- 
ronments (i.e. in voids, sheets or filaments of the Cosmic Web). 
We compare Alcyoneus to the other 151 GRGs with known lobe 
pressures in the top panel of Figure [10]'?] Alcyoneus reaffirms 
the negative correlation between length and lobe pressure 
ozy achals[002|Machaski Juro 200) and has the 

owest lobe pressures found thus far. Alcyoneus’ lobe pressure 
is in fact so low, that it is comparable to the pressure in dense 
and hot parts of the WHIM: for a baryonic matter (BM) density 


H The pysynch code is publicly available online: https://github. 
com/mhardcastle/pysynch 


We have included all publications that provide pressures, energy 
densities or magnetic field strengths. Note that some authors assume 


Ymin = 1, we assume Ymin = 10 and |Malarecki et al.|(2013) assume 


Ymin = 10°. If possible, angular lengths were updated using the LoTSS 
DR2 at 6" and redshift estimates were updated using the SDSS DR12. 


All projected proper lengths have been recalculated using our 
Collaboration et al. cosmology. When authors provided pres- 
sures for both lobes, we have taken the average. 
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Fig. 10: Of all GRGs with known lobe pressures, Alcyoneus 
is the most plausible candidate for pressure equilibrium 
with the WHIM. In the top panel, we explore the relation be- 
tween length and lobe pressure for Alcyoneus and 151 literature 
GRGs. In the bottom panel, we compare the lobe pressure of Al- 
cyoneus (green line) with the lobe pressures of the largest four 
similarly analysed GRGs (grey lines) and with WHIM pressures 
(red lines). 


pwm = 10 pcoQpmo and Twym = 10’ K, Pwyy = 410716 Pa. 
Here, pco is today’s critical density, so that o; oO gw o is today's 
mean baryon density. See the bottom panel of Figure [10] for 
a more extensive comparison between Pi, (green line) and 
Pwum (red lines). For comparison, we also show the lobe pres- 
sures of the four other thus-analysed GRGs with lp > 3 Mpc 
(grey lines). These are J1420-0545 of lp = 4.9 Mpc (Machal- 


B isi 3C 236 of lp = 4.7 Mpc (Schoenmakers et al. 
2000), J0331-7710 of lp = 3.4 Mpc (Malarecki et al.|2013) and 
B2147+816 of lp = 3.1 Mpc (Schoenmakers et al.[2000). 
Although proposed as probes of WHIM thermodynamics for 
decades, the bottom panel of Figure[10|demonstrates that even 


the largest non-cluster literature GRGs are unlikely to be in 
pressure equilibrium with their environment. Relying on results 


from the Overwhelmingly Large Simulations (OWLS) (Schaye 
fet al.[2010), Malarecki et al.| (2013) point out that baryon den- 


sities Ppm > 50 pco Qzpmo, which are necessary for pressure 
equilibrium in these GRGs (see the intersection of grey and red 
lines in the bottom panel of Figure |10}, occur in only 1% of 
the WHIM’s volume. By contrast, Alcyoneus can be in pres- 
sure equilibrium with the WHIM at baryon densities pgm ~ 
20 pco Qpmo, and thus represents the most promising inter- 
galactic barometer of its kind yet["] 


13 At Alcyoneus’ redshift, this density amounts to a baryon overden- 
sity of ~10. 


Martijn S.S.L. Oei et al.: The discovery of a radio galaxy of at least 5 Mpc 


Why do most, if not all, observed non-cluster GRGs have over- 
pressured lobes? The top panel of Figure[10|suggests that GRGs 
must grow to several Mpc to approach WHIM pressures in their 
lobes, and such GRGs are rare. However, the primary reason is 
the limited surface brightness sensitivity of all past and current 
surveys. Alcyoneus' lobes are visible in the LoTSS, but not in 
the NRAO VLA Sky Survey (NVSS) or in 
the Westerbork Northern Sky Survey (WENSS) (Rengelink et al. 
En Their pressure approaches that of the bulk of the WHIM 
within an order of magnitude. Lobes with even lower pressure 
must be less luminous or more voluminous, and thus will have 
even lower surface brightness. It is therefore probable that most 
GRG lobes that are in true pressure equilibrium with the WHIM 
still lie hidden in the radio sky. 


4. Conclusion 


1. We reprocess the LoTSS DR2, the latest version of the LO- 
FAR’s Northern Sky survey at 144 MHz, by subtracting an- 
gularly compact sources and imaging at 60” and 90” reso- 
lution. The resulting images allow us to 
explore a new sensitivity regime for radio galaxy lobes, and 
thus represent promising data to search for unknown GRGs 
of large angular length. We present a sample in forthcoming 
work. 

2. We discover the first 5 Mpc GRG, which we dub Alcyoneus. 
The projected proper length is lp = 4.99 + 0.04 Mpc, while 
the true proper length is at least Imin = 5.04 + 0.05 Mpc. 
We confidently associate the 20.8’ + 0.15’ radio structure 
to an elliptical galaxy with a jet-mode AGN detected in the 
DESI Legacy Imaging Surveys DR9: the SDSS DR12 source 
J081421.68+522410.0 at J2000 right ascension 123.590372’, 
God anon 52.402795° and spectroscopic redshift 0.24674+ 
6-107. 

3. Alcyoneus has a total luminosity density at v = 144 MHz 
of L, = 8 = 1-10? W Hz |, which is typical for GRGs (per- 
centile 45 + 3%). Alcyoneus’ host has a fairly low stellar 
mass and SMBH mass compared with other GRG hosts (per- 
centiles 25+9% and 23x 1196). This implies that — within the 
GRG population — no strong positive correlation between 
radio galaxy length and (instantaneous) low-frequency ra- 
dio power, stellar mass or SMBH mass can exist. 

4. The surrounding sky as imaged by the LoTSS, DESI Legacy 
Imaging Surveys, RASS and PSZ suggests that Alcyoneus 
does not inhabit a galaxy cluster. According to an SDSS-III 
cluster catalogue, the nearest cluster occurs at a comoving 
distance of 11 Mpc. A local galaxy number density count 
suggests that Alcyoneus instead inhabits a filament of the 
Cosmic Web. A low-density environment therefore remains 
a possible explanation for Alcyoneus' formidable size. 

5. We develop a new Bayesian model that parametrises in 
three dimensions a pair of arbitrarily oriented, optically 
thin, doubly truncated conical radio galaxy lobes with con- 
stant monochromatic emission coefficient. We then gen- 
erate the corresponding specific intensity function, taking 
into account cosmic expansion, and compare it to data as- 
suming Gaussian image noise. We use Metropolis-Hastings 
Markov chain Monte Carlo to optimise the parameters, 
and thus determine northern and southern lobe volumes of 
1.5 € 0.2 Mpc? and 1.0 € 0.2 Mpc’, respectively. In total, the 
lobes have an internal energy of ~10°> J, expelled from the 
host galaxy over a Gyr-scale period. The lobe pressures are 
4.8 + 0.3 - 10716 Pa and 4.9 + 0.6 - 107! Pa, respectively; 


these are the lowest measured in radio galaxies yet. Nev- 
ertheless, the lobe pressures still exceed a large range of 
plausible WHIM pressures. Most likely, the lobes are still 
expanding — and Alcyoneus' struggle for supremacy of the 
Cosmos continues. 
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Fig. A.1: Alcyoneus' projected proper length just exceeds 
that of J1420-0545. The probability that Alcyoneus (green) has 
a larger projected proper length than J1420-0545 (grey) 
ski et al. is 99.9%. For both GRGs, we take into account 
uncertainty in angular length and spectroscopic redshift, as well 
as the possibility of peculiar motion along the line of sight. 


Appendix A: J1420-0545 comparison 
We verify that Alcyoneus is the largest known radio galaxy (RG) 


in projection by comparing it with J1420-0545 
[2008), the literature's record holder. 

The angular lengths of Alcyoneus and J1420-0545 are @ = 
20.8’ + 0.15’ and à = 17.4’ + 0.05’, respectively. For J1420- 
0545, we adopt the angular length reported by 
because it lies outside the LoTSS DR2 coverage. The spec- 
troscopic redshifts of Alcyoneus and J1420-0545 are Zspec = 
0.24674 + 6 - 10? and Zspec = 0.3067 + 5 - 1074, respectively. 
For both giants, we assume the peculiar velocity along the line 
of sight uy to be a Gaussian random variable (RV) with mean 
0 and standard deviation 100 km s^!, similar to conditions in 
low-mass galaxy clusters. 

Equations describe how to calculate the cosmological red- 
shift RV z via the peculiar velocity redshift RV zp: 


p’ c’ p 1 - f, , 1 4 z, . 


Here, c is the speed of light in vacuo. Finally, we calculate the 
projected proper length RV J, = rg (z, Wi) - $9. Here, ry is the 
angular diameter distance RV, which depends on cosmological 
model parameters Wi. Propagating the uncertainties in angular 
length ¢, spectroscopic redshift Zspec and peculiar velocity 
along the line of sight up through Monte Carlo simulation, 
the projected proper lengths of Alcyoneus and J1420-0545 are 
ly = 4.99 + 0.04 Mpc and lp = 4.87 + 0.02 Mpc, respectively. 
We show the two projected proper length distributions in 
Figure The probability that Alcyoneus has the largest 
projected proper length is 99.9%. This result is insensitive to 
plausible changes in cosmological parameters; for example, 
the high-Hp (i.e. Ho > 70 km s! Mpc™') cosmology with M = 
(h = 0.7020, Ogy o = 0.0455, Quo = 0.2720, Q4 o = 0.7280) 
yields a probability of 99.8%. 


(A.1) 


Appendix B: Inclination angle comparison 


Under what conditions is Alcyoneus not only the largest GRG 
in the plane of the sky, but also in three dimensions? To an- 
swer this question, we compare Alcyoneus to the five previ- 
ously known GRGs with projected proper lengths above 4 Mpc, 
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Fig. A.2: When is Alcyoneus not only the largest GRG in 
the plane of the sky, but also in three dimensions? Alcy- 
oneus' inclination angle 0 is not well determined, and there- 
fore the full range of possibilities is shown on the horizontal 
axis. To surpass Alcyoneus in true proper length, a challenger 
must have an inclination angle (vertical axis) of at most Alcy- 
oneus' (grey dotted line). More specifically, as a function of 6, 
we show the maximum inclination angle for which challengers 
with a projected proper length /5;; > 4 Mpc trump Alcyoneus 
(coloured curves). The shaded areas of parameter space repre- 
sent regimes with a particularly straightforward interpretation. 
One can imagine populating the graph with five points (located 
along the same vertical line), representing the ground-truth in- 
clination angles of Alcyoneus and its five challengers. If any of 
these points fall in the red-shaded area, Alcyoneus is not the 
largest GRG in 3D. If all points fall in the green-shaded area, 
Alcyoneus is the largest GRG in 3D. 


which we dub challengers. A challenger surpasses Alcyoneus in 
true proper length when 


liz l bie.. 
P^.» —L, or sinb, < sina, (B.1) 


sins sin’ I 


Il, >I, or 


where le, |; and 6, are the challenger's true proper length, 
projected proper length and inclination angle, respectively. Be- 
cause the arcsine is a monotonically increasing function, a chal- 
lenger surpasses Alcyoneus if its inclination angle obeys 


lp. 
Os < Ogax c (0) , where Omax, (0) ‘= arcsin (3 sin o) (B.2) 


P 


In Figure [A.2] we show ga (0) for the five challengers with 
Ine € {4.11 Mpc, 4.35 Mpc, 4.60 Mpc, 4.72 Mpc, 4.87 Mpc} 
(coloured curves). Alcyoneus is least likely to be the longest 
GRG in 3D when its true proper length equals its projected 
proper length; i.e. when 0 = 90°. The challengers then surpass 
Alcyoneus in true proper length when their inclination angles 
are less than 55°, 61°, 67°, 71° and 77°, respectively. For 6 < 90°, 
the conditions are more stringent. 
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The third and fourth largest challengers, whose respec- 
tive SDSS DR12 host names are J100601.73+345410.5 and 
J093139.03+320400.1, harbour quasars in their host galaxies. If 
small inclination angles distinguish quasars from non-quasar 


AGN, as proposed by the unification model (e.g. |Hardcastle & 
ETET a, these two challengers may well be the longest 
radio ga 


axies in three dimensions. 


Appendix C: Lobe volumes with truncated double 
cone model 


Appendix C.1: Synopsis 


We build a Metropolis-Hastings Markov chain Monte Carlo 
(MH MCMC) model, similar in spirit to the model of 
for galaxy cluster halos, in order to formalise the 
determination of RG lobe volumes from a radio image. To this 
end, we introduce a parametrisation of a pair of 3D radio galaxy 
lobes, and explore the corresponding parameter space via the 
Metropolis algorithm|'*|For each parameter tuple encountered 
during exploration, we first calculate the monochromatic emis- 
sion coefficient (MEC) function of the lobes on a uniform 3D 
grid representing a proper (rather than comoving) cubical vol- 
ume. The RG is assumed to be far enough from the observer 
that the conversion to a 2D image through ray tracing simpli- 
fies to summing up the cube’s voxels along one dimension, and 
applying a cosmological attenuation factor. This factor depends 
on the galaxy’s cosmological redshift, which is a hyperparame- 
ter. We blur the model image to the resolution of the observed 
image, which is also a hyperparameter. Next, we calculate the 
likelihood that the observed image is a noisy version of the 
proposed model image. The imaged sky region is divided into 
patches with a solid angle equal to the PSF solid angle; the noise 
per patch is then assumed to be an independent Gaussian RV. 
These RVs have zero mean and share the same variance, which 
is another hyperparameter — typically obtained from the ob- 
served image. We choose a uniform prior over the full phys- 
ically realisable part of parameter space. The resulting poste- 
rior, which contains both geometric and radiative parameters, 
allows one to calculate probability distributions for many inter- 
esting quantities, such as the RG’s lobe volumes and inclination 
angle. The inferences depend weakly on cosmological parame- 
ters W. Furthermore, their reliability depends significantly on 
the validity of the model assumptions. 


Appendix C.2: Model 
Appendix C.2.1: Geometry 


We model each lobe in 3D with a truncated right circular cone 
with apex O € R?, central axis unit vector à € S? and opening 
angle y € [0, 2], as in Figure |9| The lobes share the same O, 
which is the RG host location. Each central axis unit vector can 
be parametrised through a position angle y € [0, 27) and an in- 
clination angle 0 € [0, z]. Each cone is truncated twice, through 
planes that intersect the cone perpendicularly to its central axis. 
Thus, each truncation is parametrised by the distance from the 
apex to the point where the plane intersects the central axis. 
The two inner (d;;,d;? € R29) and two outer (do, do? € Rso) 
truncation distances are parameters that we allow to vary inde- 
pendently, with the only constraint that each inner truncation 


^ The more general Metropolis—Hastings variant need not be consid- 
ered, as we work with a symmetric proposal distribution. 


distance cannot exceed the corresponding outer truncation dis- 
tance. 


Appendix C.2.2: Radiative processes 


The radiative formulation of our model is among the simplest 
possible. The radio emission from the lobes is synchrotron ra- 
diation. We consider the lobes to be perfectly optically thin: 
we neglect synchrotron self-absorption. The proper MEC is as- 
sumed spatially constant throughout a lobe, though possibly 
different among lobes; this leads to parameters jy1, jy? € Rso. 
The relationship between the specific intensity /, (in direction 
? at central frequency ve) and the MEC j, (in direction f at cos- 
mological redshift z and rest-frame frequency v = ve (1 + z)) is 


HONO 


hee MOD ® x) 
(1 +2) 


(+z) 


where / represents proper length. The approximation is valid 
for a lobe with a spatially constant MEC that is small enough 
to assume a constant redshift for it. A/(?) is the proper length 
of the line of sight through the lobe in direction 7. The inferred 
MECs j,4 (v) , jy2 (v) thus correspond to rest-frame frequency v. 


Appendix C.3: Proposal distribution 


In order to explore the posterior distribution on the parameter 
space, we follow the Metropolis algorithm. The Metropolis al- 
gorithm assumes a symmetric proposal distribution. 


Appendix C.3.1: Radio galaxy axis direction 


To propose a new RG axis direction given the current one whilst 
satisfying the symmetry assumption, we perform a trick. We 
populate the unit sphere with N € Ns, points (interpreted 
as directions) drawn from a uniform distribution. Of these N 
directions, the proposed axis direction is taken to be the one 
closest to the current axis direction (in the great-circle distance 
sense). Note that this approach evidently satisfies the criterion 
that proposing the new direction given the old one is equally 
likely as proposing the old direction given the new one. Also 
note that the distribution of the angular distance between cur- 
rent and proposed axis directions is determined solely by N. 

In the following paragraphs, we first review how to perform 
uniform sampling of the unit two-sphere. More explicitly than 
in [Scott & Tout] (1989), we then derive the distribution of the 
angular distance between a reference point and the nearest of 
N uniformly drawn other points. The result is a continuous uni- 
variate distribution with a single parameter N and finite support 
(0, 7). Finally, we present the mode, median and maximum like- 
lihood estimator of N. As far as we know, these properties are 
new to the literature. 


Uniform sampling of S? Let us place a number of points 
uniformly on the celestial sphere S?. The spherical coordinates 
of such points are given by the RVs (b, ©), where ® denotes 
position angle and O denotes inclination angle. As all posi- 
tion angles are equally likely, the distribution of ® is uniform: 
o ~ U[O0, 27). In order to effect a uniform number density, the 
probability that a point lies within a rectangle of width dy and 
height dé in the (y, 0)-plane equals the ratio of the solid angle of 
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the corresponding sky patch and the sphere’s total solid angle: 


sin 6 dy dé 


Py<®<y+dy,0<0<6+d6) = 1 
7t 


(C.2) 


The probability that the inclination angle is found somewhere 
in the interval [0, 0 + d0), regardless of the position angle, is 
therefore 


P(0 < © < 0 + d0) = dFe(0) = fo(0)d0 


Qm as 
sin 0 dO 1 
Í ae p 5 sin 


(C.3) 


where Fo is the cumulative distribution function (CDF) of 0, 
and fọ the associated probability density function (PDF). So, 


— cos Q 


1 P 1 
fe(0) = 2 sin@; Fe(0):— | fo(0^) do’ = 5 (C.4) 
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Nearest-neighbour angular distance distribution Pick a 
reference point and stochastically introduce N other points in 
above fashion, which we dub its neighbours. We now derive the 
PDF of the angular distance to the nearest neighbour (NNAD). 
Let (ger, ref) be the coordinates of the reference point and let 
(p, 0) be the coordinates of one of the neighbours. Without 
loss of generality, due to spherical symmetry, we can choose to 
place the reference point in the direction towards the observer: 
bref = 0. (Note that Yer is meaningless in this case.) The angular 
distance between two points on S? is given by the great-circle 
distance £. For our choice of reference point, we immediately 
see that E(Qrer, Oret: 9, 0) = 0. Because @ is a realisation of O, £ 
too can be regarded as a realisation of an RV, which we call =. 
Evidently, the PDF fz(£) = fo(é) and the CDF Fz(£) = Fe(£). 
Now consider the generation of N points, whose angular dis- 
tances to the reference point are the RVs (Ej) := (ZEj,..., Ey). 
The NNAD RV M is the minimum of this set: M := min{&;}. 
What are the CDF Fy and PDF fy of M? 


Fy(p) := P(M < u) = P(minimum of (Zi) < u) 
= P(at least one of the set (Ej) < y) 
= ] - P(none of the set {€;} € yu) 


= | - P(all of the set (£;) > y). (C.5) 


Because the (5j) are independent and identically distributed, 


N 
Fu(u) 21- | [PE > 


i=] 


=-1-P%(E>p)=1-(1- Fz)". (C.6) 


By substitution, the application of a trigonometric identity and 
differentiation to u, we obtain the CDF and PDF of M: 


Fu (u) = 1 — cos?" (5) fuu) = N sin (E) cogi (5). (C.7) 


In Figure|C.1| we show this PDF for various values of N. 

The mode of M (i.e. the most probable NNAD), Hmoae, is the solu- 
tion to “i (mode) = 0. The median of M, Umedian, is the solution 
to F m(Hmedian) = j. Hence, 


1 
Hmode = arccos f = x [median = arccos (27 = 1). (C.8) 


N 
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As common sense dictates, both equal a for N = 1 and tend to 
0 as N — ov. We find the mean of M through integration by 
parts: 


E[M] = Í iuis | PEE, 
0 


= [utu] - | Fu) dg 


= [ cos?" (5) du = nr cos?" (w) du. (C.9) 
0 2 0 
Again via integration by parts, 
S 2k-l n [2N 
su -»| | a = zy) (C.10) 


Maximum likelihood estimation A typical application is 
the estimation of N in the PDF fulu | N) Ln us- 
ing data. Let us assume we have measured k NNADs, denoted 
by (ji, .... Uk}. Let the joint PDF or likelihood be 


k 
LN) := | | fu (ui | N) 
i-i 
N\A 
= (=x) B sin uj (COS uj + jt 
i=l 


(C.11) 


To find Nur, we look for the value of N that maximises L(V). 
To simplify the algebra, we could however equally well max- 
imise a k-th of the natural logarithm of the likelihood, or the 
average log-likelihood Î := k^! In £(N), because the logarithm is 
a monotonically increasing function: 


^ 1 
(N) := z PLW) =InN-NIn2 


k 
1 
+ — 2 In sin u; + (N — 1) In(cos u; + 1). (C.12) 
k i=] 
We find Myr g by solving d (Nure) = 0. This leads to 
E 7 
NMLE = [na e 2; In(cos u; + 1) (C.13) 


An easy limit to evaluate is the case when 41, ..., 4g — 0. In 
such case, cos u; — 1, and so 1 pr In(cos u; + 1) > 1n2. Then, 
Nme > (04)! — oo. This is expected behaviour: when all 
measured NNADs approach 0, the number of points distributed 
on the sphere must be approaching infinity. 


Appendix C.3.2: Other parameters 


The other proposal parameters are each drawn from indepen- 
dent normal distributions centred around the current parameter 
values. These proposal distributions are evidently symmetric, 
but have support over the full real line, so that forbidden pa- 
rameter values can in principle be proposed. As a remedy, we set 
the prior probability density of the proposed parameter set to 0 
when the proposed opening angle is negative or exceeds 5 rad, 
at least one of the proposed MECS is negative, or when at least 
one of the proposed inner truncation distances is negative or 
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Fig. C.1: Probability density functions (PDFs) of the nearest-neighbour angular distance (NNAD) RV M between some fixed point 
and N other points distributed randomly over the celestial sphere. As the sphere gets more densely packed, the probability of 
finding a small M increases. For each N, we provide the mean point number density n. 


exceeds the corresponding proposed outer truncation distance. 
In such cases, the posterior probability density is 0 too, as it is 
proportional to the prior probability density. Consequently, the 
Metropolis acceptance probability vanishes and the proposal is 
rejected. We do not enter forbidden regions of parameter space. 
The condition of detailed balance is still respected: probability 
densities for transitioning towards the forbidden region are 0, 
just as probability densities for being in the forbidden region. 


Appendix C.4: Likelihood 


We assume the likelihood to be Gaussian. To avoid dimension- 
ality errors, we multiply the likelihood by a constant before we 
take the logarithm: 


Np 
in(£- (o 2x) -- A 2, (ro [i] — Tym Fi]? . (C.14) 
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Here, o is the image noise, N, € R59 is the number of reso- 
lution elements in the image, Nj € N is the number of pixels 
in the image, and /,, [i] and J,» [i] are the i-th pixel values of 
the observed and modelled image, respectively. For simplicity, 
one may multiply the likelihood by a constant factor (or, equiva- 
lently, add a constant term to the log-likelihood): the acceptance 
ratio will remain the same, and the MH MCMC runs correctly. 


Appendix C.5: Results for Alcyoneus 


We apply the Bayesian model to the 90” LoTSS DR2 image of Al- 
cyoneus, shown in the top panel of Figure|9| Thus, the hyperpa- 
rameters are z = 0.24674, ve = 144 MHz (so that v = 180 MHz), 
ÜpwuM = 90", N = 750 and o = V2. 1.16 Jy deg ?. We set 
the image noise to V2 times the true image noise to account 
for model incompleteness. This factor follows by assuming that 
the inability of the model to produce the true lobe morphol- 
ogy yields (Gaussian) errors comparable to the image noise. To 
speed up inference, we downsample the image of 2,048 by 2,048 
pixels by a factor 16 along each dimension. We run our MH 
MCMC for 10,000 steps, and discard the first 1,500 steps due to 
burn-in. 


Table [C1] lists the obtained maximum a posteriori probabil- 
ity (MAP) estimates and posterior mean and standard deviation 
(SD) of the parameters. 


Table C.1: Maximum a posteriori probability (MAP) estimates 
and posterior mean and standard deviation (SD) of the param- 
eters from the Bayesian, doubly truncated, conical radio galaxy 
lobe model of Section [3.9] 


parameter MAP estimate posterior mean and SD 
91 307° 307 21° 
p2 140° 139 + 2° 
|01 — 90°| 54° 5172 
|02 — 90°| 25° 18+7° 
yı 9° 101° 
y2 24° 26 42° 
dii 2.7 Mpc 2.6 + 0.2 Mpc 
di 4.3 Mpc 4.0 x 0.2 Mpc 
di» 1.6 Mpc 1.5 x 0.1 Mpc 
do» 2.0 Mpc 2.0 + 0.1 Mpc 
jaw) 17 Jy deg? Mpc™! | 17 + 2 Jy deg? Mpc”! 
jaQ) |22Jydeg? Mpc™! | 18 3 Jy deg? Mpc"! 


The proper volumes V; and V; are derived quantities: 


E m 3 p 
V= jum y (ds dẹ), (C.15) 
just like the flux densities F,1 (ve) and F,2 (ve) at central fre- 
quency ve: 


jy () V 


Fy v) = ——L——. 
= a9 iG 


(C.16) 


Together, V and F,(v;) imply a lobe pressure P and a magnetic 
field strength B, which are additional derived quantities that we 
calculate through pysynch. 

Table|C.2|lists the obtained MAP estimates and posterior mean 
and SD of the derived quantities. 
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Fig. C.2: Our Bayesian model yields strongly correlated es- 
timates for j, (v) and V that reproduce the observed lobe 
flux densities. We show MECs j, (v) at v = 180 MHz and 
proper volumes V of Metropolis-Hastings Markov chain Monte 
Carlo samples for the northern lobe (purple dots) and south- 
ern lobe (orange dots). The curves represent all combinations 
(jy (v) , V) that correspond to a particular flux density at the 
LoTSS central wavelength ve = 144 MHz. We show the ob- 
served northern lobe flux density (purple curve) and the ob- 
served southern lobe flux density (orange curve). 


Table C.2: Maximum a posteriori probability (MAP) estimates 
and posterior mean and standard deviation (SD) of derived 
quantities from the Bayesian, doubly truncated, conical radio 
galaxy lobe model of Section B.9] 


derived quantity | MAP estimate | posterior mean and SD 


Ay 167° 168 + 2° 
Vi 1.5 Mpc? 1.5 + 0.2 Mpc? 
Vo 0.8 Mpc? 1.0 + 0.2 Mpc? 
Fy e) 63 mJy 63 +4 mJy 
Fy5 (ve) 44 mJy 45 + 5 mJy 
Pan 4.7 - 10716 Pa 4.8 + 0.3 - 10716 Pa 
Pona 5.4 - 10716 Pa 5.0 + 0.6 - 10776 Pa 
Pa 4.8 - 10716 Pa 4.9 + 0.3 - 10716 Pa 
Paso 54.10 6Pa | 5.0+0.6-107'° Pa 
Brin, 45 pr 45 +1 pT 
Bmin,2 48 pT 46+3 pT 
Beg! 42 pr 43 +1 pT 
Beg2 45 pT 43 +3 pT 
Eini 6.3 - 10° J 6.2 + 0.4 - 10° J 
Emin2 3.7.10? J 4.4 + 0.6- 10? J 
Esa 64-10? J 6.3 + 0.4.10? J 
Eeq2 3.8.10? J 4.4 + 0.6- 10? J 


The uncertainties of the parameters and derived quantities 
reported in Tables{C-TJand(C Zare not necessarily independent. 
To demonstrate this, we present MECs and volumes from the 
MH MCMC samples in Figure [C.2] MECS and volumes do not 
vary independently, because their product is proportional to 
flux density (see Equation |C.16); only realistic flux densities 
correspond to high-likelihood model images. 


Finally, we explore a simpler variation of the model, in 
which we force the lobes to be coaxial. In such a case, the true 
proper length / and projected proper length /, are additional 
derived quantities: 


doi + do 
p= 1 7. | = [sin. 


(C.17) 
COS y 
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For Alcyoneus, this simpler model does not provide a good fit 
to the data. 


